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Abstract We analyze the light curves of the recent solar eclipses measured 
by the Herzberg channel (200-220 nm) of the Large Yield RAdiometer (LYRA) 
onboard PROBA-2. The measurements allow us to accurately retrieve the center- 
to-limb variations (CLV) of the solar brightness. The formation height of the 
radiation depends on the observing angle so the examination of the CLV provide 
information about a broad range of heights in the solar atmosphere. We employ 
the ID NLTE radiative transfer COde for Solar Irradiance (COSI) to model the 
measured light curves and corresponding CLV dependencies. The modeling is 
used to test and constrain the existing ID models of the solar atmosphere, e.g. 
the temperature structure of the photosphere and the treatment of the pseudo- 
continuum opacities in the Herzberg continuum range. We show that COSI can 
accurately reproduce not only the irradiance from the entire solar disk, but also 
the measured CLV. It hence can be used as a reliable tool for modeling the 
variability of the spectral solar irradiance. 



1. Introduction. 

The variability of the solar irradiance may have a direct impact on climate 
(see e.g. the recent reviews by |Haigh| 120071 and |Gray et aL|l20lH| . Although the 
measurements and modeling of the solar irradiance were under the close attention 
during the last decade, the complete picture of the solar variability is still far 
from being clear (see e.g. Harder et aU 120091 Haigh et aZ-t 12010]) . Therefore the 



launch of every new space mission devoted to the measurements of the solar 
irradiance is able to provide a crucial piece of complementary information as 
well as to nourish the theoretical models. 
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Figure 1. Relative variations of the irradiance as measured by the Hcrzbcrg channel of LYRA. 
The intervals of zero intensity occur during the occultations when PROBA-2 passes the Earth 
shadow in the winter season. The panels show two eclipses separated by the three occultations 
on January 15, 2010 (a) and four eclipses on July 11, 2010 (b). The periodic abrupt changes 
of the irradiance level are due to the spacecraft maneuvers. 



In this paper we analyze the first measurements of the Large Yield RA- 
diometer (LYRA) ( jHochedez et al., 2006 ; B cnmoussa et al., 2009[ ) onboard the 
PROBA-2 satellite launched on November 2, 2009. Up to now LYRA has ob- 
served several solar eclipses (see Fig. [TJ). 

During the eclipse the Moon consecutively covers different parts of the solar 
disk. The light curve of the eclipse depends on the CLV of the solar brightness 
and on the geometry of the eclipse (the angular radii of the Sun and the Moon as 
well as the minimum distance between their centers which is reached during the 
maximum phase of the eclipse). If the geometry of the eclipse is known and the 
distribution of solar brightness has radial symmetry then the light curve of the 
eclipse can be used to retrieve the CLV of the solar brightness. Let us notice that 
the assumption of the radial symmetry is well-justified for the January 15, 2010 
eclipse as the solar activity level was very low (according to the USAF/NOAA 
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data the total sunspot area was about 0.025% of the full solar disk). The CLV 
of the solar brightness provide a valuable information about the solar atmosphere 
(Allcndc Prieto, Asplund, and Fabiani Bendicho[ 1200 4; Kocste rke, Allende Prieto, and Lambert[ | 



2008} and determine the irradiance variations on the time-scale of the solar rota- 



tion ( Fligge, Solanki, and Unruh, 2000 1 . Additionally the changes of the spectral 
solar irradiance during the eclipses are important for studying the Earth's atmo- 
sphere response (K ishore Kumar, Su brahma nyam, and John}l2011IISumod et aU [ 

wni 

|Shapiro et~aL\ I2TJTU1 showed that the ID NLTE radiative transfer COdc for 
Solar Irradiance (COSI) (see |Haberreiter, Schmutz, and Huben y. 2008) allows to 
calculate the solar spectrum (125 nm-1 /im) from the entire solar disc with a high 
accuracy. In this paper we use COSI to calculate the CLV of the solar brightness 
and compare them with ones deduced from the eclipse light curves as observed 
by LYRA. We show that the measured CLV provide an important constrains 
on the UV opacities and the temperature structure of the solar atmosphere. We 
come up with a model which allows to accurately reproduce the measurements. 

We restrict ourselves to the modeling of the eclipse profiles and CLV in 
the Hcrzberg channel of LYRA. The solar irradiance in the Hcrzbcrg contin- 
uum range (200-220 nm) is of especial importance for the climate modeling 
as it directly affects the ozone concentration and stratospheric temperature 
flBrasseur et al, 1987| |Rozanov et al, 2 006 Sh apiro et al, 20iTc| ). The proper 
modeling of the formation of this radiation in the solar atmosphere is a base 
for the variability modeling and for the irradiance reconstruction to the past 
(see e.g. |Krivova, Solanki, and Unruh"} 120111 IVieira et all 120111 |Shapiro et al.\ 

We are aware that the ID models do not necessarily reflect the average physi- 
cal properties of the inhomogeneous solar atmosphere ( |Uitenbroek and Criscuoli, 2011[ ) J 
Thus the main goal of this paper is not to learn the new facts about the dynamic 
3D solar atmosphere but rather to develop a reliable semi-empirical tool for 
modeling the solar irradiance variability. 

In Sect. [H we deduce the empirical CLV from the LYRA observations of the 
January 15, 2010 eclipse. In Sect. [3] we compare these CLV with ones calcu- 
lated by COSI and discuss the constrains on the temperature structure of the 
solar atmosphere (Sect. l3~2]) and UV opacities (Sect. l3~3|) . The main results are 
summarized in Sect. [4] 



2. Empirical center-to-limb variations as deduced from the LYRA 
observations 

PROBA-2 evolves on a dawn-dusk heliosynchronous orbit, with an altitude of 
720 km on average, which allows a quasi-permanent observation of the Sun. The 
spacecraft accomplishes a full orbit in about 100 minutes. In case of solar eclipse, 
it is therefore not unusual that it crosses the eclipse zone more than once. Lyra 
data in Hcrzbcrg channel are acquired, for all three units, by experimental PIN 
detectors made of diamond. Such detectors have proven to be very stable with re- 
spect to temperature variations ( |Benmoussa et al, 2004[ ). The nominal cadence 
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Figure 2. The profiles of the 05:00 (a) and 09:00 (b) transits retrieved from the level 1 data 
as well as the deviations between these profiles and the profile retrieved from the level 2 data 
of the unit 2 (c,d). The shaded area on panel c indicates the estimated error range. The dotted 
lines correspond to the widest part of the error range in the time interval between 5 h and 
5.3 h (i.e. excluding the problematic feature after 5.3 h). 



of acquisition is 20 Hz. A more detailed discussion of the in- flight performance 
of LYRA is given in |Dominique et al] 120111 

The first light curve of the eclipse event was obtained by LYRA on January 
15, 2010. The eclipse was shortly preceded by LYRA first light on January 6, 
2010 and was the longest annular solar eclipse of the millennium. It was observed 
on the ground from Africa and Asia and was seen as a partial from the PROBA- 
2. The eclipse lasted more than 6 hours, so the PROBA-2 passed through 
the Moon's shade three times. However the intermediate transit could not be 
observed due to the simultaneous occultation (i.e. it was shaded by the Earth). 

The raw (level 1) data collected by the Herzberg channel of LYRA during 
this eclipse are presented in panels a and b of Fig. [2] The plotted data were 
corrected for the dark current which is still present in the original data. The 
05:00 UTC transit of the January 15, 2010 eclipse (hereafter first transit) was 
simultaneously observed by the LYRA units 1 and 2 (the back-up and standard 
acquisition units, accordingly), while the 09:00 UTC transit (hereafter second 
transit) only by the unit 2. The drops of the irradiance after the first and before 
the second transits correspond to the occultations (see panel a of Fig. [T]). 
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The level 1 data are uncalibrated. The calibrated (level 2) data are also 
available for the community. For the analysis presented below, we will therefore 
use the profiles retrieved from the level 2 data. These data are corrected for the 
temperature effects, degradation and the dark current. Let us note that level 
2 data always refer to the unit 2 measurements, while the measurements from 
the back-up units 1 and 3 are normally used for the calibration and currently 
only available in their uncalibrated version ( |Dominique et al., 2011[ ). Originally 
the level 2 data were corrected for the degradation by adding a time-dependent 
offset (to allow a better analysis of the solar flares). The offset shifts the zero- 
level of the irradiance and leads to the erroneous profiles of the eclipses and 
solar variability ( |Shapiro et al., 2011d I . Therefore for our analysis the offset was 



removed from the level 2 data. 

The blue curves on the panels c and d of Fig. [2] represent the deviations 
between the level 1 and level 2 data. The zero- level of both datasets was corrected 
as discussed above. One can see that the level 2 data are almost identical to the 
original measurements of the standard acquisition unit 2. This confirms that the 
temperature correction for the diamond detectors is very small. From now on 
we will not distinguish between the level 1 and level 2 data of the unit 2. On the 
contrary there is a significant deviation between the data collected by the unit 
1 and unit 2 during the first transit. Both units were carefully calibrated in the 
ground facilities so their different responses are connected with the degradation 
which affects the sensitivities of the units. The standard acquisition unit 2 was 
opened almost constantly, while the back-up units 1 and 3 were opened only 
occasionally. As a result on January 15, 2010 unit 2 degraded approximately by 
15% while unit 1 did not yet show any noticeable degradation. The differential 
behavior of the units allows us to estimate the error range of the measured 
profiles. We defined the uncertainty that it is twice the difference between the 
profiles as observed by the unit 2 and the unit 1 (see panels c of Fig. [2]) . This 
approach does not include any systematic deviations, common to both units but 
it measures the degradation due to the space exposure. The arbitrary character 
of this estimate does not affect the results presented below. The second transit 
was only observed by the unit 2 so for its analysis we will use the error range 
defined from the first transit. 

The profile of the irradiance variation during the eclipse depends on the 
CLV of the solar brightness. For example a strong decrease of the brightness 
towards the solar limb would lead to a smaller residual irradiance during the 
maximal phase of the eclipse and accordingly to a deeper eclipse profile. Thus, 
the observed profiles allow to determine the CLV. For simplicity we adopted the 
widely used polynomial parametrization of the CLV (see e.g. Ncckcl and Labs. 
HMNeckelll2nn5) : 

a * m 
7 = 2^ A ^ , i 1 ) 

-'center 

% 

where /center is the disk-center intensity, [i is the cosine of the heliocentric angle 
and Ai are free parameters. The solar disk was divided into the thirteen suppos- 
edly uniform concentric rings and the brightness of each ring was calculated with 
the help of the Eq. ((TJ using the \i value of the ring's mean circle. We checked 
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Figure 3. (a) Empirical CLV deduced from the profile of the first transit as observed by 
the unit 2 (orange curve); (b) the deviations between the CLV plotted on panel (a) and the 
CLV deduced from the profile of the first transit as observed by the unit f (red curve) and 
CLV deduced from the profile of the second transit as observed by the unit 2 (blue curve). The 
empirical CLV deduced from the two extreme profiles corresponding to the edges of the error 
range in Fig. [2] delimit the CLV error range (the shaded area). 



that such division of the solar disc allows to calculate the irradiance profiles with 
the accuracy better than 0.02% and thus is sufficient for our purposes. Simple 
geometrical calculations allow to obtain the profile of the eclipse for each set 
of the coefficients, assuming that the time dependencies of the angular distance 
between the Sun and the Moon as well as of their angular sizes are known. 

For the both transits of the January 15, 2010 eclipse we searched for the set 
of the coefficients Ai which lead to the best agreement with the observed profiles 
and minimized the error £ defined as: 

^ ^ ] (^obs ^emp) s 

i 

where P bs is the observed profile and "P C mp is the empirical profile calculated 
from the Eq. (JTJ) . For the first and second transits the summation was done over 
the time intervals between 5 h and 5.3 h and between 9 h and 9.25 h respectively. 

The minimization was performed applying the additional condition of the 
monotonous CLV. The CLV profiles were calculated using the second degree 
polynomial. The employment of the higher degree polynomial (up to sixth de- 
gree) had no visible effect. The CLV depend on the wavelength ( |Neckel, 1 996; 
Hcstroffcr and Magnan, 1998). Thus, the coefficients Ai determined by the min- 
imization procedure correspond to the CLV of the solar intensity convolved with 
the profile of the LYRA Herzberg channel. The latter is a combined profile of 
the detector and filter. 

The resulting empirical CLV dependencies are presented in Fig. [3] The CLV 
dependency for the first transit is normalized to unity in the disk center, while 
all other CLV dependencies are normalized to give the same integral flux from 
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Figure 4. The deviations between the profiles as measured by unit 2 and as calculated with 
the empirical CLV. The error bars (dotted lines) are the same as in Fig. [2] 



the entire disk. The error range of the CLV was estimated performing the mini- 
mization procedure to the two maximal error profiles of the first transit. So the 
shaded area in the Fig. [3] is constrained by the two CLV dependencies which 
correspond to the edges of the shaded area in the panel c of Fig. [5] 

The deviations between the profiles as measured by LYRA and calculated 
with the empirical CLV are shown in Fig. 2J The deviations can be attributed 
to the limited accuracy of the measurements (note that they are within the 
estimated error bars) and to the violation of the radial symmetry of the solar 
brightness. One of the sources of the asymmetry of the solar brightness is the 
inhomcgenous structure of the quiet Sun, which consists from several bright- 
ness components (Fontenl a et ai, 1999| ). The small amplitude of the deviations 
supports the analysis presented below. 



3. Comparison with modeling 

To calculate the theoretical CLV we employed the ID NLTE radiative transfer 
code COSI developed by |Hubeny[ 119811 IHamann and Schmutzl 119871 



Schmutz, Hamann, and Wessolowski [T989 Habcrrci ter, Schmutz, and Hubeny] 



2008] |Shapiro et al.\\201Q[ COSI simultaneously solves the statistical equilibrium 
and radiative transfer equations in the spherically symmetrical geometry. The 
temperature and density structures of the different components of the solar 
atmosphere were taken from IFontenla et alX 119991 while the electron density 
and all level populations were self-consistently calculated in the NLTE. 

3.1. Solar irradiance in the Herzberg continuum range 

The proper calculation of the radiative transfer in the Herzberg continuum region 
is sophisticated by two factors. Firstly, the continuum opacity in this region is 



strongly affected by the NLTE overionization (see e.g. Shchukina and Trujillo Bueno[ [ 
120011 IShort and H auschildtl I2005P , so the proper NLTE calculations are neces- 
sary. Secondly, the immense number of weak, mostly spectrally unresolved lines, 
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form the so-called UV line haze in this region. The creation of the correct line list 
which would include all possible spectral lines is a task of a tremendous difficulty. 
Although considerable progress was reached during the last few decades, there 
is still only 1% of the UV lines which are measured in the laboratory, while all 
remaining lines are predicted only theoretically flKurucz, 2005[ ). As the conse- 
quence the existing line list are not complete and underestimate the opacity in 
the UV spectral region (Busa et al, 200l[|Haberreiter, Schmutz, and Hubeny, 2008|| 



Short and Hauschildt, 2009[ ). 

To account for the missing opacity |Shapiro et al. [[2"0 101 multiplied the contin- 
uum opacity coefficient fc c (A) by the wavelength dependent coefficient / C (A): 

k' c (X) = k c (X) • / C (A). (2) 

The employed multiplicative coefficient is a step function of the wavelength with 
the step equals to 1 nm. It was empirically determined so that the UV irradiance 
calculated by COSI equals to the irradiance as measured by SOLSTICE (SOLar- 
STellar Irradiance Comparison Experiment, see McClintock, Snow, and Woods[ 



I2TJCJ5)) onboard the SORCE satellite ( |Rottman, 2005[ ) during the 2008 solar mini 



mum. The additional opacity which is necessary to reproduce the SOLSTICE/SORCEJ 
irradiance does not exceed a few percents of the total opacity included in COSI. 
It is only necessary in the 160-320 nm spectral region. Shortward of 160 nm the 
total opacity is mainly dominated by the photo-ionization continuum opacity, 
while longward of 320 nm the existing line lists are accurate enough to reproduce 
the irradiance with better than 5% accuracy (and with better than 2% accuracy 
longward of 400 nm). 

As during the solar minimum the solar irradiance (except the extreme UV) is 
dominated by the quiet Sun, the solar atmosphere model C (average supergranule 
cell interior model) from IFontenla et at, 1999|was used for the calculations. Let 
us however note that it would be possible to be in agreement with the SOL- 
STICE/SORCE measurements while using another temperature profile by ad- 
justing additional opacities. Thus the SOLSTICE/SORCE measurements alone 
docs not allow to constrain the temperature structure of the solar atmosphere. 
In Sect. I3.2l we show that although simultaneous adjustment of the temperature 
structure and additional opacities does not affect the irradiance from the entire 
solar disc, it significantly changes the CLV. The comparison with the measured 
CLV allows to choose the most suitable for the CLV calculations model. 

|Shapiro et aL|[2T)10l assumed that the coefficient of the additional opacity does 
not depend on height in the solar atmosphere. This implies that it docs not have 
any dominant source as the concentration of every particular ion or molecule 
is height dependent. In Sect. 13.31 we show that the CLV provide important 
complementary information which allow to reevaluate this assumption and help 
to understand the nature of the additional opacity. 

3.2. Test of the temperature structure 

To test the sensitivity of the CLV to the change of the temperature structure 
we performed the calculations employing the models for several different compo- 
nents of the solar atmosphere: model A (faint supergranule cell interior), model 
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C (average supergranule cell interior), model E (quiet network) and model P 
(plage). The temperature and density structures were taken from lFontenla et al.[ 
1999, The models A and E correspond to the cold and warm components of 
the quiet Sun, while the Model C represents the spatially averaged quiet Sun. 
Sha piro et aL|l2"0"T0l showed that the calculations with the latter model can repro- 
duce spectral irradiance measured by SOLSTICE (up to 320 nm) and SIM (Solar 
Irradiancc Monitor; Harder et al\ I2005|) (from 320 nm onward) onboard the 
SORCE satellite during the 2008 solar minimum, as well as SOLSPEC (SOLar 
SPECtral Irradiance Measurements; IThuillier et ail I2004p during the ATLAS 
3 mission in 1994 with high accuracy. They used the model C to calculate the 
/c(A) factor (see Eq. (J2J)) for the additional opacities in the 160-320 nm spectral 
range. 

The same procedure of the / C (A) factor fitting was performed for the Mod- 
els A and E. The Model A is colder than Model C so it yields smaller UV 
irradiance. Thus a smaller additional opacity is necessary to reproduce the SOL- 
STICE/SORCE measurements in the 160-320 nm spectral region. Accordingly 
the use of Model E leads to a larger additional opacity. The Model P yields 
so high irradiance that it is not possible to reach SOLSTICE/SORCE level by 
increasing the continuum opacity. Thus we didn't recalculate the / C (A) factor 
for the model P and left it the same as for the Model C. 

In Fig. [S]we present the calculated CLV for each of these models. The radi- 
ation which comes from the regions close to the solar limb is formed in higher 
and colder regions of the solar atmosphere than the radiation coming from the 
disc center. Hence the solar brightness is decreasing towards the limb. One can 
see that the colder the model, the stronger CLV it yields. This can be partly 
explained by the fact that the sensitivity of the Planck function to the tem- 
perature change is the decreasing function of the temperature. Thus the same 
change of temperature results in larger alteration of the Planck function and 
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Table 1. Theoretical profiles vs. observed and empir- 
ical (calculated with the empirical CLV) profiles. The 
numbers are discrepancies calculated with the help of 
Eq. J3]l. The minimum values of the discrepancy are 
boldfaced. 
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49.9 
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216.3 


209.2 


245.2 


236.3 



accordingly larger change of the emergent irradiance for the colder models. The 
complete picture depends on the temperature and density structure as well as on 
the opacity behavior in each of the solar atmosphere components. It is further 
sophisticated by the NLTE effects which causes deviations of the source function 
from the Planck function. 

Although due to the readjustment of the additional opacity the calculations 
with the Models A and E yield the same UV irradiance as the calculation with 
the Model C, the corresponding CLV dependencies are remarkably different. All 
differences between the theoretical and empirical CLV's have a sudden drop at 
R/Rsun ~ 0.99 (see panel (b) of Fig. [5]). The rings which correspond to these 
points (see Sect. [2]) have very small relative area, so the eclipse profiles are 
basically insensitive to the change of their brightness. Thus the reliability of the 
empirical CLV for these points is very low. 

Interestingly, the calculations with the model C underestimate the CLV and 
are outside of the estimated error region. One of the possible reasons for this 
could be the erroneous assumption of the depth independency of the additional 
opacity coefficient f c (see Eq. ©). This possibility will be discussed in the 
Sect. 13.31 On the other hand it is known that the ID models can underestimate 



the anisotropy of the radiation field ( Shapiro et al., 2011b| Kleint et al., 2011) 
and the CLV of the solar brightness ( |Koesterke, Allende Prieto, and Lambert, 20"08} | 



Uitenbroek and Criscuoli, 2011 ). Hence the reported disagreement could be a 
signature of the general problems which are inherent to ID modeling. 

The warmer Models E and P yield even weaker CLV than the Model C. At 
the same time the colder model A yields the CLV which are in a good agreement 
with the empirical ones. 

The corresponding deviations between the calculated and measured profiles 
are shown in Fig. [5] The structure common to the both transits represents 
properties of the models, while individual structure is due to the limited accuracy 
of the measurements and violation of the radial symmetry of the solar brightness. 
We can define the discrepancy between the profiles V\ and V\ as: 

A 1 , 2 = W- 4 .^J2(Vi-n) 2 /Nj, (3) 
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Figure 6. The same as Fig.[4]but with the eclipse's profiles calculated for the solar atmosphere 
models A, C, E and P. 



where N is number of the time points in the profiles and the scaling factor 10 -4 
is introduced for convenience. 

In Table [T] we present the discrepancies between the theoretical, measured 
and empirical (calculated with the empirical CLV, see Sect. [5]) profiles. One can 
see that the calculations with the Model A arc in the good agreement with the 
measurements and is very close to the empirical profile. Let us however note 
that the model A is not able to properly reproduce the visible and near infrared 
irradiance as well as the main molecular bands (e.g. CH G-band and CN violet 
system) in the solar spectrum. The calculations with the Models E and P result 
in too weak CLV and consequently in too high residual irradiance during the 
maximal phase of the eclipse. Although the calculations with the Model C are 
outside of the estimated error range the yielded profiles are reasonably close to 
the observations (e.g. they will be hardly distinguishable from the measured ones 
in the scale of the Fig.[T]). Thus the Model C still can be used in the calculations 
when the high accuracy is not necessary. 

3.3. Adjustments of the additional opacities 

The introduction of the additional opacity described in Sect. 13.11 changes the 
formation heights of the UV radiation and thus affects the CLV. All calculations 
presented in Sect. 13.21 were performed assuming depth independent coefficient f c 
in Eq. ^ . At the same time if the additional opacity arises from the unaccounted 
lines of some particular atom or molecule X then it should be scaled with the 
relative concentration of X. Then the coefficient f c becomes height dependent: 



f c (\,h) = l + T c (\) 



nx(h)/n tota .i(h) 
{n x {h)/n tota x{h)) 1 



(4) 



where h is the height in the solar atmosphere, n to tai and nx are the total and 
species X concentrations respectively. J- c {\) is a step function which can be 
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determined empirically in the same way as the f c factor was determined in 
Sha piro et al.\ 120101 The height dependency of the additional opacity is intro- 
duced by the fraction in the right side of the equation. The latter is the relative 
concentration of the species X, normalized for convenience to unity. 

To better understand the possible origin of the additional opacity we scaled 
it with the relative concentrations of CN and CO molecules and Fe I ion. In 
each case we empirically determined the T C (X) factor to reproduce the SOL- 
STICE/SORCE UV measurements. In panel a of Fig. [7] we show the dependency 
of the fraction from the Eq. @ on height. The CN and CO molecules are 
mainly present in the narrow photospheric layer. As the CO molecule has larger 
dissociation potential than the CN (11.1 eV against 7.76 eV) it is more sensitive 
to the temperature change ( |Berdyugina, Solanki, and Frutiger, 2003 ) and the 



maximum of the CO relative concentration is very close to the temperature 
minimum. The maximum of the CN relative concentration is slightly shifted to- 
wards the lower levels of the photosphere where the effect of the density increase 
overcompensate the temperature increase. Iron is mostly ionized throughout the 
solar atmosphere. Similarly to the molecular case the neutral iron concentration 
has a peak slightly below the temperature minimum. At the same time it does 
not drop so abruptly in the lower level of the solar photosphere and starts to 
increase in the chromosphere due to the strong density decrease. 

On the three lower panels of Fig. [7] we show the contribution functions 
for the intensities in the Herzberg continuum range. They were obtained by 
convolving the individual contribution functions (see, e.g. |Gray| 119921 p. 151) 
for every frequency with the combined profile of the detector and filter in the 
Herzberg channel. Thus the plotted contribution functions sample the regions 
where the irradiance as measured by the Herzberg channel of LYRA is formed. 
The scaling of the additional opacity with the Fe I or molecular concentrations 
alters the optical depth scale and simultaneously shifts the limb and disc center 
contribution functions to the higher and lower levels of the solar atmosphere 
respectively. This increases the distance between the peaks of the limb and disc 
center contribution functions, which results in the stronger CLV. The effect is 
the most prominent for the scaling with the CO relative concentration (due to 
the strong shift of the limb contribution function) and the least prominent for 
the Fe I case (due to a small effect on the disc center contribution function). 

The above discussion is confirmed by the Figs. |S] and M where the CLV depen- 
dencies and the eclipse profiles for the different cases of the additional opacity 
scaling are presented. One can see that the scaling of the additional opacity 
makes the CLV stronger. The scaling with the CO relative concentration results 
in too strong CLV and accordingly too small residual irradiance during the 
maximal phase of the eclipse. At the same time the scaling with the Fe I or CN 
relative concentrations moves the CLV dependencies and irradiance profiles very 
close to the observed ones and hence can solve the problem addressed in the end 
of Sect. 13.21 Let us notice that, in opposite to the calculations with the Model 
A, these calculations can also properly reproduce all main features in the solar 
spectrum. 

While the relative variations of the irradiance can be measured by LYRA 
with a very high precision, the absolute calibration of LYRA radiometers is 
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Figure 7. Opacity scaling factor (the fraction in the right side of the Eq. J4{) and temperature 
(black curve) as a function of height (a) as well as the contribution functions for near the limb 
(R/Rsun = 0.95) and disc center intensities. The contribution functions are shown for the 
case of height independent additional opacity (black curves in (b), (c) and (d)) and additional 
opacity scaled with the CN (b), CO (c) and Fe I (d) relative concentrations. The zero-point 
of the height scale is denned as the layer at which the continuum optical depth at 500 nm is 
equal to one. 
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Figure 8. The same as Fig. [3] but for different modes of calculations with the Model C. 
The "CN" , "CO" and "Fe I" points were calculated scaling the additional opacity with the 
relative concentration of the corresponding species. The "0.8" and "1.5" points were calculated 
adjusting the additional opacity so that COSI yields 80% and 150% of the UV irradiancc 
measured by the SOLSTICE/SORCE. The "C" points correspond to the standard calculations 
using Model C (i.e. with height independent additional opacities adjusted to reproduce the 
SOLSTICE/SORCE UV measurements). 



Table 2. The same as Table [TJ but for different modes of cal- 
culations with the Model C and for the altered levels of the 
SOLSTCIE/SORCE measurements. The minimum values of 
the discrepancy for different modes of calculations with the 
Model C arc boldfaced. 





Passa 


ge 1 


Passaj 


;e 2 




observed 


empirical 


observed 


empirical 


Standard 


105.0 


90.9 


129.3 


111.5 


CN 


47.8 


19.1 


66.5 


25.0 


CO 


83.4 


73.3 


85.4 


60.8 


Fc I 


59.5 


35.0 


83.9 


55.0 


80% 


170.8 


161.7 


198.8 


187.4 


150% 


51.4 


15.6 


64.0 


3.1 



very tricky and renders the determination of the absolute level of the irradi- 
ance almost impossible without the use of an external reference. The additional 
opacity in the COSI were adjusted to reproduce the irradiance as measured by 
SOLSTCIE/SORCE. If the solar UV irradiance were different from the mea- 
sured by the SOLSTCIE/SORCE, then the different additional opacity would 
be necessary to reproduce it. The readjustment of the additional opacities will 
affect not only the absolute level of the measured irradiance but also the CLV. 
Thus the CLV deduced from the eclipse analysis can be used to indirectly test 
the SOLSTCIE/SORCE measurements. With this goal we made two experi- 
ments readjusting the additional opacity so that COSI yields 80% and 150% of 
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the SOLSTCIE/SORCE irradiance (the same scaling factor was applied for all 
wavelengths in the 160-320 nm interval). In these experiments we followed the 
approach of |Shapiro et q/4l2010l and assumed that the f c factor from the Eq. ^ 
is independent on height. 

The corresponding CLV are given on Fig. [5J while the eclipse profiles are 
plotted on the bottom panels of Fig. [5] One can see that the increase of the 
absolute level of the irradiance leads to a stronger CLV and good agreement 
with the measurements. Let us note that due to the problem of the standard 
modeling discussed in Sect. 13.21 this should be considered not as a contradiction 
to the SOLSTICE/SORCE measurements but rather as a consequence of too 
weak CLV yielded by the calculations with the Model C and depth independent 
coefficient f c . At the same time the decrease of the irradiance level leads to a 
significantly weaker CLV and strong deviations from the measurements. 

In Table [5] we compare profiles which were discussed above with the LYRA 
measurements. The first line corresponds to the calculations with height in- 
dependent additional opacity adjusted to reproduce the SOLSTICE/SORCE 
measurements. 
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4. Conclusions 

The profiles of the eclipse light curves provide an important information for 
testing and refining the solar atmosphere models. We have shown that the eclipse 
profiles observed by the Herzberg channel of LYRA are in a very good agreement 
with the synthetic profiles calculated with the ID NLTE radiative transfer code 
COSI. 

The calculated profiles are very sensitive to the temperature structure of the 
solar atmosphere and to the treatment of the UV opacity, whose significant 
part can still be missing from the modern models. The best agreement between 
the observed and measured profiles can be reached in two different regimes of 
calculations. The first one corresponds to the calculations with the Model A 
atmosphere structure from IFontenla et ail 119991 and constant coefficient of the 
additional opacity. The second one is the calculations with the Model C but 
assuming that the missing opacity scales with the relative concentration of the 
neutral iron or CN molecule. Let us notice that while the calculations with the 
Model A are not able to properly reproduce the near UV, visible and infrared 
solar irradiance, the calculations with the model C yield a good agreement with 
measurements over the entire solar spectrum ( |Shapiro et al., 2010[ ). 

Our results could hint that the missing opacity originates in the layer a few 
hundred kilometers below the temperature minimum and could be due to the 
unaccounted lines of the neutral iron or another element with the similar ioniza- 
tion potential (e.g. silicon or magnesium) or due to the unaccounted molecular 
lines (e.g. CN). 

We are aware of the limitations of modeling with ID solar atmosphere. 
|Koesterke, Allende Prieto, and Lambert) l2008l and lUitenbroek and Criscuolii 1201 1 | 
show that 3D modeling can lead to a stronger CLV of the solar brightness. 
This prevents us from making an unambiguous conclusion. Nevertheless the ID 
radiative transfer code are widely used for the interpretation of the the stellar and 
solar spectra and presently they are the de facto standards for modeling the solar 
spectral irradiance variability ( |Krivova and Solanki, 2008 Fontenla et at, 2 009; 
IDomingo et al., 20090 . 

The variability of the solar irradiance on the 11-year and solar rotational 
time-scales is usually attributed to the competition between the irradiance in- 
crease due to the bright active components (e.g. plage and active network) and 
irradiance decrease due to the dark sunspots (see e.g. lKrivova et aO 120030 . The 
contrast between active features and the quiet Sun strongly depends on the disk 



position (Unruh, Solanki, and Fligge, 1999 Fligge, Solanki, and Unruh, 2000). 



Thus the center-to-limb variations of the solar brightness analyzed in this paper 
play an important role in the modeling of the solar irradiance variability. The fact 
that the calculations with COSI arc in a good agreement with the measurements 
strongly supports its suitability for such modeling ( |Shapiro et al., 2011d| ). 

Acknowledgements The research leading to this paper was supported by the Swiss Na- 
tional Science Foundation under grant CRSI122-130642 (FUPSOL) and grant 200020-130102. 
We thank the LYRA PROBA 2 science team for their work in producing the data sets used 
in this paper and their helpful recommendations. 



SOLA: shapiroetal_eclipses.tex; 1 February 2012; 4:03; p. 16 



Eclipses observed by LYRA 



References 

Allende Prieto, C, Asplund, M., Fabiani Bendicho, P.: 2004, Ccnter-to-limb variation of solar 
line profiles as a test of NLTE line formation calculations. Astron. Astrophys. 423, 1109 — 
1117. doi: 10.1051/0004-6361:20047050. 

Benmoussa, A., Schiihle, U., Haenen, K., Nesladek, M., Koizumi, S., Hochedez, J.-F.: 2004, PIN 
diamond detector development for LYRA, the solar VUV radiometer on board PROBA II. 
Physica Status Solidi Applied Research 201, 2536-2541. doi: 10. 1002/pssa. 200405187 

Benmoussa, A., Dammasch, I.E., Hochedez, J.-F., Schiihle, U., Roller, S., Stockman, Y., 
Scholzc, F., Richter, M., Kroth, U., Laubis, C, Dominique, M., Krctzschmar, M., Mekaoui, 
S., Gissot, S., Theissen, A., Giordanengo, B., Bolsee, D., Hermans, C, Gillotay, D., Dcfisc, 
J.-M., Schmutz, W.: 2009, Pre-fiight calibration of LYR A, the solar VUV radiometer on 
board PROBA2. Astron. Astrophys. 508, 1085-1094. doi 10.1051/0004-6361/200913089 

Bcrdyugina, S..V., Solanki, S..K., Frutiger, C: 2003, The molecular zeeman effect and diagnos- 
tics of solar and stellar magnetic fields, ii. synthetic stokes profiles in the zeeman regime. 
Astron. Astrophys. 412, 513. 

Brasscur, G., de Rudder, A., Keating, G.M., Pitts, M.C.: 1987, Response of Middle Atmosphere 
to Short-Term Solar Ultraviol et Variations. II. Theory. J. Geophys. Res. 92, 903-914. 
doi:10.1029/JD092iD01p00903 

Busa, I., Andretta, V., Gomez, M.T., Terranegra, L.: 2001, A method to estimate the effect of 
line blanketing in NLTE radiative transfer calculations. Astron. Astrophys. 373, 993-997. 
doi: 10. 1051/0004-6361:20010661 

Domingo, V., Ermolli, 1., Fox, P., Frohlich, C, Haberrcitcr, M., Krivova, N., Kopp, G., 
Schmutz, W., Solanki, S.K., Spruit, H.C., Unruh, Y., Vogler, A.: 2009, Solar Surface 
Magnetism and Irradiancc on Time Scales from Days to the 11- Year Cycle. Space Sci. 
Rev. 145, 337-380. doi: 10.1007/sll214-009-9562-l 

Dominique, M., Hochedez, J.-F., Schmutz, W., Dammasch, I., BenMoussa, A., Shapiro, A. I., 
Krctzschmar, M.: 2011, The LYRA instrument on-board PROBA2: description and in-flight 
performances. Solar Phys. This issue.. 

Fligge, M., Solanki, S.K., Unruh, Y.C.: 2000, Modelling irradiance variations from the surface 
distribution of the solar magnetic field. Astron. Astrophys. 353, 380-388. 

Fontcnla, J.M., Curdt, W., Haberreiter, M., Harder, J., Tian, H.: 2009, Semiempirical Models of 
the Solar Atmosphere. III. Set of Non-LTE Models for Far-Ultraviolct/Extremc-Ultraviolct 
Irradiance Computation. Astrophys. J. 707, 482-502. doi|T0~1088/0004-637X/ 707/l/482| 

Fontenla, J., White, O.R., Fox, P.A., Avrett, E.H., Kurucz, R.L.: 1999, Calculation of 
Solar Irradianccs. I. Synthesis of the Solar Spectrum. Astrophys. J. 518, 480-499. 
doi: 10.1086/307258 

Gray, D.F.: 1992, The observation and analysis of stellar photospheres. 

Gray, L.J., Beer, J., Geller, M., Haigh, J.D., Lockwood, M., Matthes, K., Cubasch, U., Fleit- 
mann, D., Harrison, G., Hood, L., Lutcrbacher, J., Meehl, G.A., Shindell, D., van Geel, 
B., White, W.: 2010, SOLAR INFLU ENCES ON CLIMATE. Reviews of Geophysics 48, 
RG4001. doi:10.1029/2009RG000282 

Haberreiter, M., Schmutz, W., Hubcny, I.: 2008, NLTE model calculations for the solar atmo- 
sphere with an iterative treatment of opacity distribution functions. Astron. Astrophys. 
492, 833-840. doi [17110 5 1 /0004-636 1 : 2008095*031 

Haigh, J.D.: 2007, The Sun and the Earth's Climate. Living Reviews in Solar Physics 4, 2-. 

Haigh, J.D., Winning, A.R., Toumi, R., Harder, J.W.: 2010, An influence of solar spectral 
variations on radiative forcing of climate. Nature 467, 696-699. doi 10.1038/nature09426 



Hamann, W.-R., Schmutz, W.: 1987, Computed He II spectra for Wolf-Rayet stars - A grid of 

models. Astron. Astrophys. 174, 173-182. 
Harder, J.W., Fontenla, J., Lawrence, G., Woods, T., Rottman, G: 2005, The Spectral Ir- 
radiancc Monitor: Measurement Equations and Calibration. Solar Phys. 230, 169-204. 

doi:10.1007/sll207-005-1528-l 
Harder, J.W., Fontenla, J.M., Pilewskie, P., Richard, E.C., Woods, T.N.: 2009, Trends in solar 

spectral irradiancc variability in the visible and infrared. Geophys. Res. Lett. 36, L07801. 

doi:10.1029/2008GL036797 
Hcstroffcr, D., Magnan, C: 1998, Wavelength dependency of the Solar limb darkening. Astron. 

Astrophys. 333, 338-342. 
Hochedez, J.-F., Schmutz, W., Stockman, Y., Schiihle, U., Benmoussa, A., Koller, S., Haenen, 

K., Berghmans, D., Defise, J.-M., Halain, J. -P., Theissen, A., Dclouillc, V., Slcmzin, V., 




SOLA: shapiroetal_eclipses.tex; 1 February 2012; 4:03; p. 17 



Shapiro et al. 



Gillotay, D., Fussen, D., Dominique, M., Vanhcllemont, F., McMullin, D., Kretzschmar, 
M., Mitrofanov, A., Nicula, B., Wautcrs, L., Roth, H., Rozanov, E., Riicdi, I., Wehrli, C, 
Soltani, A., Amano, H., van der Linden, R., Zhukov, A., Clette, F., Koizumi, S., Mortct, 
V., Remes, Z., Petersen, R., Nesladek, M., D'Olieslaeger, M., Roggen, J., Rochus, P.: 
2006, LYRA, a solar UV radiometer on Proba2. Advances in Space Research 37, 303 — 312. 
doi: 10. 1016/j.asr.2005. 10.041 
Hubeny, I.: 1981, Non-LTE analysis of the ultraviolet spectrum of A type stars. II Theoretical 
considerations and interpretation of the VEGA Lyman-alpha region. Astron. Astrophys. 
98, 96-111. 

Kishore Kumar, K., Subrahmanyam, K.V., John, S.R.: 2011, New insights into the strato- 
spheric and mesosphere-lower thermospheric ozone response to the abrupt changes in solar 
forcing. Annates Geophysicae 29, 1093-1099. doi 10.5194/angeo-29-1093-2011 

Klcint, L., Shapiro, A.I., Berdyugina, S.V., Bianda, M.: 2011, Solar turbulent magnetic fields: 
Non-LTE modeling of the Hanle effect in the C2 molecule. Astron. Astrophys.. submitted. 

Koesterke, L., Allende Prieto, C, Lambert, D.L.: 2008, Center-to-Limb Variation of So- 
lar Three-dimensional Hydrodynamical Simulations. Astrophys. J. 680 , 764 — 773. 
doi: 10.1086/587471 

Krivova, N.A., Solanki, S.K.: 2008, Models of solar irradiancc variations: Current status. 

Journal of Astrophysics and Astronomy 29, 151-158. doi 10.1007/sl2036-008-0018-x 
Krivova, N.A., Solanki, S.K., Unruh, Y.C.: 2011, Towards a long-term record of solar total and 

spectral irradiancc. Journal of Atmospheric and Solar- Terrestrial Physics 73, 223-234. 

doi:10.1016/j.jastp.2009.11.013 
Krivova, N.A., Solanki, S.K., Fliggc, M., Unruh, Y.C.: 2003, Reconstruction of solar irradiancc 

variations in cycle 23: Is solar surface magnetism the cause? Astron. Astrophys. 399, 1-4. 

doi j| 10. 1051/0004-636 1 : 20030029] 
Kurucz, R.L.: 2005, Including all the lines. Memorie delta Societa Astronomica Italiana 

Supplement 8, 86 — . 

McClintock, W.E., Snow, M., Woods, T.N.: 2005, Solar-Stellar Irradiance Comparison Ex- 
periment II ( SOLSTICE II): Pre-Launch and On-Orbit Calibrations. Solar Phys. 230, 
259-294. doi |10.10077sli 207-005-1585-5 

Ncckcl, H.: 1996, On the wavelength dependency of solar limb darkening (AA303 to 1099 nm). 
Solar Phys. 167, 9-23. 

Neckel, H.: 2005, Analytical Reference Functions F(A) for the Sun's Limb Darkening and 
Its Absolute Continuum Inte nsities (AA 300 to 1100 m). Solar Phys. 229, 13-33. 
doi: 10.1007/sll207-005-4081-z 

Neckel, H., Labs, D.: 1994, Solar limb darkening 1986-1990 (AA303 to 1099nm). Solar Phys. 
153, 91. 

Rottm an, C: 2005, The SORCE Mission. Solar Phys. 230, 7-25. 
doi:10.1007/sll207-005-8112-6 

Rozanov, E., Egorova, T., Schmutz, W., Peter, T.: 2006, Simulation of the stratospheric 
ozone and temperature response to the solar irradiancc variability during sun ro- 
tation cycle. Journal of Atmospheric and Solar- Terrestrial Physics 68, 2203-2213. 
doi:10.1016/j.jastp.2006.09.004 

Schmutz, W., Hamann, W.-R., Wcssolowski, U.: 1989, Spectral analysis of 30 Wolf-Rayet stars. 
Astron. Astrophys. 210, 236-248. 

Shapiro, A. I., Schmutz, W., Schoell, M., Haberreiter, M., Rozanov, E.: 2010, NLTE 
solar irradiance modeling with the COSI code. Astron. Astrophys. 517, A48+. 
doi |rO)51/0004-63 61/200913987l 

Shapiro, A.I., Schmutz, W., Rozanov, E., Schoell, M., Haberreiter, M., Shapiro, A.V., Nyeki, S.: 
2011a, A new approach to the long-term reconstruction of the solar irradiancc leads to large 
historical solar forcing. Astron. Astrophys. 529, A67+. doi: 10.1051/0004-6361/201016173 

Shapiro, A. I., Fluri, D.M., Berdyugina, S.V., Bianda, M., Ramelli, R.: 2011b, NLTE modeling 
of Stokes vector ccntcr-to-limb variations in the CN violet system. Astron. Astrophys. 
529, A139+. doi: 10.1051/0004-6361/200811299 

Shapiro, A.V., Rozanov, E., Egorova, T., Shapiro, A. I., Peter, T., Schmutz, W.: 2011c, Sensi- 
tivity of the Earth's middle atmosphere to short-term solar variability and its dependence 
on the choice of solar irradiancc data set. Journal of Atmospheric and Solar- Terrestrial 
Physics 73, 348-355. doi |10.1016/j.jastp.2010.02.011 

Shapiro, A.V., Shapiro, A.I., Schmutz, W., Dominique, M., Dammasch, I., Wehrli, C: 2011d, 
Solar rotational cycle as observed by LYRA. Solar Phys. This issue.. 



SOLA: shapiroetal_eclipses.tex; 1 February 2012; 4:03; p. 18 



Eclipses observed by LYRA 



Shchukina, N., Trujillo Bueno, J.: 2001, The Iron Line Formation Problem in Three- 
dimensional Hydrodynamic Models of Solar-like Photospheres. Astrophys. J. 550, 
970-990. doi: 10.1086/319789 

Short, C.I., Hauschildt, P.H.: 2005, A Non-LTE Line -Blanketed Model of a Solar-Type Star. 
Astrophys. J. 618, 926-938. doi: 10.1086/426128 

Short, C.I., Hauschildt, P.H.: 2009, Non-LTE Model ing of the Near-Ultraviolet Band of Late- 
Type Stars. Astrophys. J. 691, 1634-1647. doi 10.1088/0004-637X/691/2/1634 

Sumod, S.G., Pant, T.K., Vinceth, C, Hossain, M.M., Antonita, M.: 2011, Response of the 
tropical mesopause to the longest annular solar eclipse of this millennium. Journal of 
Geophysical Research (Space Physics) 116, A06317. doi: 10.1029/2010JA016326] 

Thuillier, G., Floyd, L., Woods, T.N., Cebula, R., Hilsenrath, E., Herse, M., Labs, D.: 2004, 
Solar Irradiancc Reference Spectra. In: J. M. Pap, P. Fox, C. Frohlich, H. S. Hudson, 
J. Kuhn, J. McCormack, G. North, W. Sprigg, & S. T. Wu (ed.) Solar Variability and its 
Effects on Climate. Geophysical Monograph 141, Washington DC American Geophysical 
Union Geophysical Monograph Series 141, 171 — . 

Uitenbroek, H., Criscuoli, S.: 2011, Why One-dimensional Models Fail in the Diagnosis of 
Average Spectra from Inhomogcneous Stellar Atmospheres. Astrophys. J. 736, 69 — . 
doi: 10. 1088/0004-637X/736/1/69 

Unruh, Y.C., Solanki, S.K., Fliggc, M.: 1999, The spectral dependence of facular contrast and 
solar irradiancc variations. Astron. Astrophys. 345,635-642. 

Vieira, L.E.A., Solanki, S.K., Krivova, N.A., Usoskin, I.: 2011, Evolution of the solar irradiancc 
during the Holocene. Astron. Astrophys. 531, A6+. doi: 10.1051/0004-6361/201015843 



SOLA: shapiroetal_eclipses.tex; 1 February 2012; 4:03; p. 19 



SOLA: shapiroetal_eclipses.tex; 1 February 2012; 4:03; p. 20 



